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In a kagome lattice, the time reversal symmetry can be broken by a staggered magnetic flux 
emerging from ferromagnetic ordering and intrinsic spin-orbit coupling, leading to several well- 
separated nontrivial Chern bands and intrinsic quantum anomalous Hall effect. Based on this idea 
and ab initio calculations, we propose the realization of the intrinsic quantum anomalous Hall effect 
in the single layer Cs 2 Mn 3 Fi 2 kagome lattice and on the (001) surface of a Cs 2 LiMn 3 Fi 2 single 
crystal by modifying the carrier coverage on it, where the band gap is around 20 meV. Moreover, 
a simplified tight binding model based on the in-plane dda antibonding states is constructed to 
understand the topological band structures of the system. 

PACS numbers: 71.20.-b, 73.43.-f, 71.70.Ej, 73.20.At 


The quantum anomalous Hall (QAH) insulator is 
known as a topological state of matter in two dimen¬ 
sions (2D) with dissipationless chiral edge states pro¬ 
tected by the electron band topology [1-19], which can 
be used to design new quantum devices such as the chi¬ 
ral interconnect [20]. Haldane proposed a QAH model 
of a 2D honeycomb lattice with a staggered magnetic 
flux that averages out to zero [1]. Recently, the QAH 
effect has been theoretically predicted in magnetic topo¬ 
logical insulators [3-6] and observed experimentally in Cr- 
doped (Bi,Sb) 2 Te 3 (CBST) thin films at a low temper¬ 
ature around 30 mK [10]. However, such a low working 
temperature severely hinders the practical application of 
the QAH effect, which is largely constrained by the spa¬ 
tial inhomogeneity induced by multiple random dopants 
in the system [21]. For practical applications, to reach 
a higher working temperature, it is important to search 
for intrinsic QAH materials without random doping and 
increase the topologically nontrivial band gap, as well as 
the ferromagnetic (FM) Curie temperature (Tc). 

As is known, a honeycomb lattice with spin-orbit cou¬ 
pling (SOC) and FM ordering naturally realizes the Hal¬ 
dane model and the intrinsic QAH effect, and there are 
a lot of proposals for the QAH effect based on mag¬ 
netic doped graphene (silicene) or graphene/FM (anti¬ 
ferromagnetic) heterostructures [22-26]. However, stoi¬ 
chiometric magnetic honeycomb compounds exist rarely 
in nature. On the other hand, the kagome lattice [Fig. 
1(a)] shares similar topological physics with the honey¬ 
comb lattice [27, 28], and there are a number of magnetic 
materials adopting a kagome lattice structure [29-31]. 
With a staggered magnetic flux emerging from SOC and 
the FM ordering, the electronic bands of a 2D kagome lat¬ 
tice exhibit nontrivial Chern numbers, making it possible 
to realize the intrinsic QAH effect in such systems. As 
an illustrative example, we propose here that the QAH 
effect can be realized in a single layer Cs 2 Mn 3 Fi 2 kagome 
lattice and on the (001) surface of a Cs 2 LiMn 3 Fi 2 single 
crystal by modifying the surface carrier concentration so 
that the surface bands are at 2/3 filling. Compared to 



FIG. 1. (Color online) (a) Schematic of the Mn kagome lattice 
and the antibonding d orbitals near the Fermi level. An effec¬ 
tive staggered magnetic flux 0 in the triangles and —20 in the 
hexagons of the kagome lattice emerge due to the FM ordering 
and the intrinsic SOC. (b) The unit cell of Cs 2 LiMn 3 Fi 2 with 
space group R3, which contains three primitive cells, (c). The 
crystal structure of single layer Cs 2 Mn 3 Fi 2 . 

CBST, such materials without any magnetic doping are 
easier to synthesize, and the much higher homogeneity 
and the large band gap around 20 meV enables a much 
higher working temperature. A simplified tight binding 
(TB) model based on the in-plane dda antibonding states 
is constructed to understand the topological physics of 
the system. While most studies of the kagome lattice are 
focused on the flat band physics at 1/3 filling [32-37], 
the system considered here is at 2/3 filling with a low 
energy theory governed by Dirac fermions at K and K' 
points similar to that in the Haldane model. This class 
of materials opens a new direction for the exploration of 
the topological states of matter. 

Experimentally, the Cs 2 LiMn 3 Fi 2 single crystal is syn¬ 
thesized within the R3 space group (point group Sq) [29], 
where the Mn atoms are surrounded by a corner-sharing 
octahedron constructed by F atoms and form perfect 
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kagome lattices layer by layer as shown in Fig. 1(b) and 
Fig. 1(a). Each [Mn 3 Fi 2 ] kagome layer is sandwiched 
by two layers of Cs atoms and connected by Li atoms 
between [Cs 2 Mn 3 Fi 2 ] layers along c axis (z direction). 
We note that the [MnFe] octahedra are strongly elon¬ 
gated by the Jahn-Teller distortion with ratios of long 
(within the layer) to short axes (out of the layer) of 
about 1.18 and form an antiferrodistortive Jahn-Teller 
ordering. This is very similar to the FM manganese flu¬ 
oride CsMnF 4 , where the antiferrodistortive ordering of 
Jahn-Teller elongated [MnFe] octahedra with the ratio 
1.17 makes the asymmetric a superexchange interactions 
between the empty orbitals and half-filled dx 2 _y 2 or¬ 
bitals of Mn^+ favor the planar FM coupling with Tc = 
19 K [38]. Because of such a large Jahn-Teller distor¬ 
tion, it is the trigonal rather than the octahedral field 
that dominates the splitting of the Sd orbital of Mn as 
shown below. 

Our ab initio calculations are carried out by the den¬ 
sity functional method [39, 40] based on the plane- 
wave ultrasoft pseudopotential approximation [41] with 
the generalized gradient approximation of the Perdew- 
Burke-Ernzerhof type for the exchange-correlation poten¬ 
tial [42], as implemented in the ESTATE (Beijing Sim¬ 
ulation Tool of Atomic Technology) package [43]. Eully 
optimized lattice constants a = 7.5698 A, c = 17.4362 
A with the force smaller than 0.01 eV/A are used in 
all calculations, which are very close to the experimen¬ 
tal measurements a = 7.440 A, c = 17.267 A. A slab 
constructed by three layers of Cs 2 LiMn 3 Ei 2 with more 
than 18 Aof vacuum are used for the surface study, and 
the virtual crystal approximation (VGA) [44, 45] is used 
to simulate the carrier concentration on the surface. The 
cutoff energy for wave function expansion is set to 400 eV, 
and 8x8x8 and 8 x 8 x 2 k meshes are used for the bulk 
and slab self-consistent calculations, respectively. SOC 
effect is considered consistently in the calculations. All 
results except VGA calculations are well reproduced by 
VASP [46, 47]. 

Our calculations conclude that the 3d states of Mn^+ 
are strongly spin polarized and Gs 2 LiMn 3 Ei 2 is a EM 
insulator with magnetic moment of 4.0 /i^/Mn. The cal¬ 
culated total energy of the EM solution is 6.71 eV/f.u. 
lower than the nonmagnetic solution and 1.57 eV/f.u. 
lower than the all-in-all-out antiferromagnetic (AEM) so¬ 
lution, which confirms that the antiferrodistortive order¬ 
ing of Jahn-Teller elongated [MnEe] octahedra favors an 
in-plane EM exchange coupling between the magnetic 
moment of Mn. This is quite different from the spin liquid 
system herbertsmithite ZnGu 3 (OH) 6 Gl 2 , where the Gu 
atoms also form the layered kagome lattice but the mag¬ 
netic moment of Gu favors an in-plane AEM exchange 
coupling [30, 48]. 

In Eig. 2(b), we show the spin resolved density of states 
(DOS) and projected density of states (PDOS) on Mn for 
the EM solution, in which the states around the Eermi 



FIG. 2 . (Color online) (a) The first Brillouin zones of 
Cs2LiMn3Fi2 with space group R 3 . The four types inversion- 
invariant points are r(0,0,0), Z( 7 r, 7 r, 7 r), F(71,71, 0 ) and 
L( 7 r, 0 , 0 ). The blue hexagon shows the 2 D Brillouin zones 
for the single layer Cs2Mn3Fi2 and slab, in which the high- 
symmetry k points f, M and K are labeled, (b) Spin re¬ 
solved DOS and PDOS of the bulk Cs2LiMn3Fi2. (c) FM 
band structures of Cs2LiMn3Fi2 without SOC, in which red 
circles, blue squares and brown diamonds denote the projec¬ 
tions on the d^2, d^2_y2^dxy, and dyz^dzx orbitals of Mn, 
respectively, (d) FM band structures of Cs2LiMn3Fi2 with 
SOC, in which the parities of the topmost three valence bands 
at eight inversion-symmetric points are marked. The inset is 
the enlargement around the K point to show the fully opened 
gap. All Fermi levels are defined at 0 eV. 


level (0 eV) are mostly contributed by Mn-3d up spin 
states, while all the down spin Mn-3(i states are pushed 
above more than 1 eV away due to the large Hand’s rule 
coupling. With the trigonal field effect taken into ac¬ 
count, the up spin Mn-3(i states roughly split into three 
groups: dz‘ 2 , dx 2 _y 2 Szdxy, and dyz&idzx^ as represented by 
the red circles, blue squares and brown diamonds in Fig. 
2(c), respectively. The dx 2 -y 2 &i,dxy group further splits 
into the bonding states and the antibonding 

states d *2 , 2 &d* , due to the in-plane Mn-Mn dda inter- 

X y J^y 

action. As a result, the three valence states and three 
conduction states are (i* 2 _^ 2 &d*^ f and d^^ respec¬ 
tively, as shown in Fig. 2(b) and Fig 2(c). In Fig. 2(c), 
there are two important features to be pointed out. (1) 
The dispersion along the 2 : direction is very weak (see the 
bands along V — Z, F — L), indicating that Cs 2 LiMn 3 Fi 2 
has a weak interlayer hopping and is therefore a good 
layer by layer material. (2) In the absence of SOC, the 
highest two valence bands cross with each other near the 
K points, and form six Weyl points protected by the crys¬ 
tal symmetry. Once the SOC is included, an insulating 
gap opens at the Weyl points, resulting in a topological 
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nontrivial phase if the Fermi level is tuned to the gray 
dashed line as shown in Fig 2(d), z.e., one hole doped per 
primitive cell. The parities of the Bloch wave functions 
of the topmost three valence bands at eight inversion- 
symmetric points r(0, 0, 0), ^(tt, tt, tt), 3F(7r, tt, 0), 3I/(7r, 
0, 0) are calculated and labeled in Fig. 2(d). Based on 
the parity analysis, if one electron per primitive cell is 
reduced from the system, the occupied bands below the 
Fermi level [the gray dashed line in Fig. 2(d)] will have 
weak topological numbers (z^o, ^ 2 , ^ 3 ) = (0, 1, 1, 1), 

which indicates the system can be viewed as many 2D 
QAH insulator layers weakly stacked along (111), i.e., 
the z direction. 

However, it is too hard to reduce one electron di¬ 
rectly from the 3D bulk of the material. Instead, it 
is much easier to modify the carrier concentration on 
the surface or in the 2D limit (single layer) for such 
kinds of layered materials [49, 50]. Experimentally, sin¬ 
gle layer Cs 2 Mn 3 Fi 2 as shown in Fig. 1(c) can be 
cleaved from the Cs 2 LiMn 3 Fi 2 single crystal. After 
the cleaving, two kinds of terminations, with or with¬ 
out the outermost Li layer, may be realized. Never¬ 
theless, since the outermost Li atoms are highly mobile 
and active, their concentration can be tuned depend¬ 
ing on the experimental conditions. For example, we 
calculated the total energy of the fluorine molecule F 2 , 
LiF, the single layer Li 2 Cs 2 Mn 3 Fi 2 and Cs 2 Mn 3 Fi 2 , and 
found that E(Cs 2 Mn 3 Fi 2 )+ 2 £'(LiF) is 1.78eV lower than 
E(Li 2 Cs 2 Mn 3 Fi 2 )+E(F 2 ), which clearly shows that the 
outermost Li atoms can be removed in the F 2 -rich envi¬ 
ronment. On the other hand, because the main function 
of the Li atoms in Cs 2 LiMn 3 Fi 2 is to provide carriers, re¬ 
moving them changes little the electronic structures and 
the stability of the lattice. To test this, we have op¬ 
timized the structure of single layer Cs 2 Mn 3 Fi 2 with a 
vacuum region thicker than 15 A, and have confirmed 
that all of its structure parameters are nearly the same 
as those of the bulk Cs 2 LiMn 3 Fi 2 . Furthermore, our 
calculations show that single layer Cs 2 Mn 3 Fi 2 also fa¬ 
vors a FM solution with a total magnetic moment 3.6 
/i^/Mn, which is 5.85 eV/(single layer) lower than the 
nonmagnetic solution and 1.58 eV/(single layer) lower 
than the all-in-all-out AFM solution. The calculated FM 
band structures of single layer Cs 2 Mn 3 Fi 2 without (blue) 
and with (red) SOC are shown in Fig. 3(b), where the 
three bands near the Fermi level are still mainly from the 
d* 2 _^ 2 &d*^ t orbitals. 

To understand the topological properties of single layer 
Cs 2 Mn 3 Fi 2 , we construct here a TB model of the kagome 
lattice based on the three antibonding states 

near the Fermi level. The corresponding antibonding ba¬ 
sis wave functions of the three Mn atoms in one unit cell 
are shown in Fig. 1(a), which are related to each other 
by ±27r/3 rotations. Meanwhile, the electron spins are 
polarized due to the FM ordering. In the presence of an 
intrinsic SOC, the nearest neighbor spin polarized TB 



FIG. 3. (Color online) (a) TB band structures with ti = 
0.25 eV, t 2 = 0 (blue) and t 2 = 0.4ti (red), respectively, 
(b) FM band structures of single layer Cs 2 Mn 3 Fi 2 with (red) 
and without (blue) SOC calculations, (c) Enlargement of (b) 
around the K point to show the Dirac point (without SOC) 
and topologically nontrivial gap (with SOC) clearly, (d) Cal¬ 
culated edge state for the semi-infinite boundary of the 2D 
kagome lattice Cs 2 Mn 3 Fi 2 . All Fermi levels are defined at 0 
eV. 


model takes the form 

n = , ( 1 ) 

{ij) 


where Qq, (cj^,) is the annihilation (creation) operator 
of an electron of spin a {a =t:i) on site i, ti, t 2 are 
real, and Uij = -hl(— 1 ) if j ^ i is hopping counter¬ 
clockwise (clockwise) about the triangle containing site 
i and site j. ti > 0 represents the direct hopping am¬ 
plitude, while t 2 comes from the intrinsic SOC term 
^{ij)a^ 5 where is the Pauli matrix for 

the electron spin [32, 51]. As a result, there is a magnetic 
flux (j) = 3 arctan(t 2 /ti) in the triangles and —20 in the 
hexagons of the kagome lattice, respectively, as shown in 
Fig. 1(a). We note that such a staggered magnetic flux 
averages out to zero and preserves the inversion symme¬ 
try. This is quite different from another class of kagome 
lattice models relying on FM ordering and Rashba SOC 
that breaks the inversion symmetry [28]. 

The physics of this TB model is quite similar to that of 
the Haldane model. The model consists of two dispersive 
upper bands |izi), \u 2 ) and a lower flat band \us). When 
^2=0 (the flux 0 = 0 ), the upper band \ ui ) and the 
middle band \u 2 ) touch linearly at the K and K' points, 
forming two Dirac cones with opposite helicities, while 
the middle band \u 2 ) and the lower flat band \us) touch 
quadratically at the F point. When ^2 > 0 (the flux 
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0 < 0 < tt), the three bands \ui)^\u 2 )^\u^) are detached 
from each other in energies, and carry Chern numbers 
—respectively. To be explicit, we have plotted 
the band structures for the TB model with ^ 2/^1 = 0 
and 12 /ti = 0.4, respectively in Fig. 3(a), where we have 
set ti = 0.25 eV. We note that the QAH effect can be 
realized at both 1/3 filling and 2/3 filling. While most 
existing theoretical works are focused on 1/3 filling fiat 
band physics and the related fractional QAH effect, it is 
hard to obtain a full insulating band gap at 1 /3 filling in 
the presence of hopping beyond nearest neighbors in real 
materials [Fig. 3(b)]. On the contrary, it is much easier 
to open an insulating band gap at 2/3 filling as shown in 
Fig. 3(c) and Fig. 4. 

For small ^ 2 , we can write down an effective Hamilto¬ 
nian in the continuous limit: 

Ti = -ifiVF'lp^CTxTzdx + 0-ydy)xp + rUsoci’^CTzTztp , ( 2 ) 

where 7 / is a four component wave function, cFx,y,z are 
Pauli matrices describing the lattice pseudospin, and 

= ±1 describes states at the K and K' points, respec¬ 
tively. The mass term rrisoc is given by rusoc = 2^2/3. 
Such a low energy Hamiltonian has exactly the same 
form as that of the Haldane model [1], which makes it 
much easier to understand the topological phase transi¬ 
tion. When rrisoc > 0, each of the two Dirac cones con¬ 
tributes a TT Berry phase flux, leading to a total Chern 
number C = +1, while, when rrisoc < 0, both of the 
Berry phase fluxes change sign and the total Chern num¬ 
ber becomes C = —1. 

As shown in Fig. 3(b) and Fig. 3(c), our calcula¬ 
tions confirm that single layer Cs 2 Mn 3 Fi 2 without SOC 
is a gapless semimetal with the Fermi level exactly cross¬ 
ing two Dirac points at K and i?', similar to that in 
graphene. When SOC is introduced into the system, an 
insulating gap of 20 meV is opened and the system be¬ 
comes a QAH insulator with Chern number C = 1, as 
demonstrated by our TB model. We have carried out the 
calculations of edge states by constructing the Green’s 
functions [52] for the semi-infinite boundary based on the 
maximally Localized Wannier function method [53, 54]. 
The results are shown in Fig. 3(d), in which one topologi¬ 
cally protected chiral edge state is present in between the 
valence band and the conduction band, clearly indicating 
the system has Chern number (7 = 1. 

Another proposal to achieve the required Fermi filling 
to realize the QAH effect in Cs 2 LiMn 3 Fi 2 is to modify the 
carrier concentration on the surface. The layered crys¬ 
tal structure guarantees that the sample can be easily 
cleaved, and the carrier concentration on the surface can 
be tuned by the coverage of the topmost Li and Cs atoms. 
By passivating the lower surface, we optimize the struc¬ 
ture of the upper surface (slab) without the topmost Li 
layer. Our calculations confirm that the structure vari¬ 
ation of the surface layer Cs 2 Mn 3 Fi 2 is negligible, and 
the hole only goes into the surface layer, keeping all the 




FIG. 4 . (Color online) (a) FM band structures of the slab 
constructed by three layers of Cs2LiMn3Fi2 with SOC. The 
inset is the enlargement around K point to show the topo¬ 
logically nontrivial gap (4 meV) clearly, (b) Calculated edge 
state for the slab constructed by three layers of Cs2LiMn3Fi2. 
The inset is the enlargement around the Fermi level to show 
the insulating gap (4 meV) and chiral edge state clearly. All 
Fermi levels are defined at 0 eV. 


other bulk layers insulating as shown in Fig. 4(a). The 
calculated magnetic moments are 3.6 /i^/Mn for the sur¬ 
face layer, and 4.0 /i^/Mn for the bulk layers. As shown 
in Fig. 4(a), there are only two bands coming from the 
^x‘ 2 -y 2 ^d%y t states of the topmost Mn layer crossing the 
Fermi level, while all the other (i *2 2 &d* , t bands of the 

^ X —y ' 

bulk Mn layers are fully occupied. Therefore, the surface 
layer of Mn behaves like a pure 2D kagome lattice, in 
which the QAH effect can be achieved with suitable car¬ 
rier concentration. The band structures for the slab with 
SOC exhibits an insulating gap of about 4 meV on the 
surface, as shown in the inset of Fig. 4(a). The calcu¬ 
lated chiral edge state of the slab is shown in Fig. 4(b), 
showing that the QAH effect with Chern number (7 = 1 
is realized on the surface of the slab. 

In conclusion, we have studied the possibility of realiz¬ 
ing the intrinsic QAH effect in materials with a ferromag¬ 
netic kagome lattice structure, and have predicted that 
Cs 2 LiMn 3 Fi 2 is such a good candidate. We constructed 
a simplified TB model with ferromagnetic ordering and 
SOC based on the in-plane dda antibonding states, and 
demonstrated the realization of the QAH effect at 2/3 
filling. Based on ab initio calculations, we showed that 
the single layer Cs 2 Mn 3 Fi 2 ferromagnetic kagome lattice 
has exactly 2/3 Fermi filling, and realizes an intrinsic 
QAH insulator with a band gap of around 20 meV. Such 
nondoped materials are much easier to synthesize and 
much more homogeneous than CBST, therefore enables 
a much higher working temperature than that in CBST. 
Further, our calculations suggested that a similar Fermi 
filling and QAH effect can also be achieved on the (001) 
surface layer of a Cs 2 LiMn 3 Fi 2 single crystal via a mod¬ 
ification of the surface carrier coverage. 

We acknowledge valuable discussions with Xiao-Liang 
Qi. This work is supported by the U.S. Department 
of Energy, Office of Basic Energy Sciences, Division 
of Materials Sciences and Engineering under Contract 
No. DE-AC02-76SE00515, by EAME, one of six centers 










5 


of STARnet, a Semiconductor Research Corporation pro¬ 
gram sponsored by MARCO and DARPA. 


[1] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[2] X. L. Qi, Y. S. Wu, and S. C. Zhang, Phys. Rev. B 74, 
045125 (2006). 

[3] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008). 

[4] C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, 
Phys. Rev. Lett. 101, 146802 (2008). 

[5] R. Li, J. Wang, X. L. Qi, and S. C. Zhang, Nature Phys. 
6, 284 (2010). 

[6] R. Yu, W. Zhang, H. J. Zhang, S. C. Zhang, X. Dai, and 
Z. Fang, Science 329, 61 (2010). 

[7] G. Xu, H. Weng, Z. Wang, X. Dai, and Z. Fang, Phys. 
Rev. Lett. 107, 186806 (2011). 

[8] D. Xiao, W. Zhu, Y. Ran, N. Nagaosa, and S. Okamoto, 
Nat. Commun. 2, 596 (2011). 

[9] A. Riiegg and G. A. Fiete, Phys. Rev. B 84, 201103 

( 2011 ). 

[10] C.-Z. Ghang, J. Zhang, X. Feng, J. Shen, Z. Zhang, 
M. Guo, K. Li, Y. On, P. Wei, L.-L. Wang, Z.-Q. Ji, 
Y. Feng, S. Ji, X. Ghen, J. Jia, X. Dai, Z. Fang, S.-G. 
Zhang, K. He, Y. Wang, L. Lu, X.-G. Ma, and Q.-K. 
Xue, Science 340, 167 (2013). 

[11] J. Wang, B. Lian, H. Zhang, Y. Xu, and S.-G. Zhang, 
Phys. Rev. Lett. Ill, 136801 (2013). 

[12] J. Wang, B. Lian, H. Zhang, and S.-G. Zhang, Phys. 
Rev. Lett. Ill, 086803 (2013). 

[13] H. Zhang, J. Wang, G. Xu, Y. Xu, and S.-G. Zhang, 
Phys. Rev. Lett. 112, 096804 (2014). 

[14] K. F. Garrity and D. Vanderbilt, Phys. Rev. B 90, 121103 
(2014). 

[15] J. Wang, B. Lian, and S.-G. Zhang, Phys. Rev. B 89, 
085106 (2014). 

[16] X. Kou, S.-T. Guo, Y. Fan, L. Pan, M. Lang, Y. Jiang, 
Q. Shao, T. Nie, K. Murata, J. Tang, Y. Wang, L. He, 
T.-K. Lee, W.-L. Lee, and K. L. Wang, Phys. Rev. Lett. 
113, 137201 (2014). 

[17] J. G. Gheckelsky, R. Yoshimi, A. Tsukazaki, K. S. 
Takahashi, Y. Kozuka, J. Faison, M. Kawasaki, and 
Y. Tokura, Nat. Phys. 10, 731 (2014). 

[18] G. Xu, J. Wang, G. Felser, X.-L. Qi, and S.-G. Zhang, 
Nano Letters 15, 2019 (2015). 

[19] G.-Z. Ghang, W. Zhao, D. Y. Kim, H. Zhang, B. A. Assaf, 
D. Heiman, S.-G. Zhang, G. Liu, M. H. W. Ghan, and 
J. S. Moodera, Nature Mater. 14, 473 (2015). 

[20] X. Zhang and S.-G. Zhang, Proc. SPIE 8373, 837309 

( 2012 ). 

[21] Y. Feng, X. Feng, Y. On, J. Wang, G. Liu, L. Zhang, 
D. Zhao, G. Jiang, S.-G. Zhang, K. He, X. Ma, Q.-K. 
Xue, and Y. Wang, Phys. Rev. Lett. 115, 126801 (2015). 

[22] Z. Qiao, S. A. Yang, W. Feng, W.-K. Tse, J. Ding, 
Y. Yao, J. Wang, and Q. Niu, Phys. Rev. B 82, 161414 
( 2010 ). 

[23] H. Zhang, G. Lazo, S. Bliigel, S. Heinze, and 
Y. Mokrousov, Phys. Rev. Lett. 108, 056802 (2012). 


[24] M. Ezawa, Phys. Rev. Lett. 109, 055502 (2012). 

[25] Z. Qiao, W. Ren, H. Ghen, L. Bellaiche, Z. Zhang, A. H. 
MacDonald, and Q. Niu, Phys. Rev. Lett. 112, 116404 
(2014). 

[26] Z. Wang, G. Tang, R. Sachs, Y. Barlas, and J. Shi, Phys. 
Rev. Lett. 114, 016603 (2015). 

[27] K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. 
B 62, R6065 (2000). 

[28] Z.-Y. Zhang, Journal of Physics: Gondensed Matter 23, 
365801 (2011). 

[29] U. Englich, G. Erommen, and W. Massa, Journal of Al¬ 
loys and Gompounds 246, 155 (1997). 

[30] M. P. Shores, E. A. Nytko, B. M. Bartlett, and D. G. 
Nocera, Journal of the American Ghemical Society 127, 
13462 (2005). 

[31] Z. E. Wang, N. Su, and E. Liu, Nano Letters 13, 2842 
(2013). 

[32] E. Tang, J.-W. Mei, and X.-G. Wen, Phys. Rev. Lett. 
106, 236802 (2011). 

[33] J. W. E. Venderbos, M. Daghofer, and J. van den Brink, 
Phys. Rev. Lett. 107, 116401 (2011). 

[34] S. Nishimoto, M. Nakamura, A. O’Brien, and P. Eulde, 
Phys. Rev. Lett. 104, 196401 (2010). 

[35] D. Green, L. Santos, and G. Ghamon, Phys. Rev. B 82, 
075104 (2010). 

[36] Q. Liu, H. Yao, and T. Ma, Phys. Rev. B 82, 045102 

( 2010 ). 

[37] Y. Nakata, T. Okada, T. Nakanishi, and M. Kitano, 
Phys. Rev. B 85, 205128 (2012). 

[38] W. Massa and M. Steiner, Journal of Solid State Ghem- 
istry 32, 137 (1980). 

[39] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 
(1964). 

[40] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

[41] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[42] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[43] Z. Pang and K. Terakura, Journal of Physics: Gondensed 
Matter 14, 3001 (2002). 

[44] L. Nordheim, Ann. Phys 9, 607 (1931). 

[45] L. Bellaiche and D. Vanderbilt, Phys. Rev. B 61, 7877 

( 2000 ). 

[46] G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 

[47] G. Kresse and J. Eurthmiiller, Phys. Rev. B 54, 11169 
(1996). 

[48] H. O. Jeschke, E. Salvat-Pujol, and R. Valenti, Phys. 
Rev. B 88, 075106 (2013). 

[49] H. Weng, G. Xu, H. Zhang, S.-G. Zhang, X. Dai, and 
Z. Pang, Phys. Rev. B 84, 060408 (2011). 

[50] M.-Y. Yao, L. Miao, N. L. Wang, J. H. Dil, M. Z. Hasan, 
D. D. Guan, G. L. Gao, G. Liu, D. Qian, and J.-f. Jia, 
Phys. Rev. B 91, 161411 (2015). 

[51] G. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

[52] M. L. Sancho, J. L. Sancho, J. L. Sancho, and J. Rubio, 
Journal of Physics E: Metal Physics 15, 851 (1985). 

[53] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997). 

[54] 1. Souza, N. Marzari, and D. Vanderbilt, Phys. Rev. B 
65, 035109 (2001). 



